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SUMMARY 


An investigation was made of the pressure and equi librium- temperature 
distributions on a sphere at Mach numbers of 2.00 and 4.15. The local 
aerodynamic heat transfer was also measured on a sphere at a Mach number 
of 2.00 and on a hemisphere -cylinder at a Mach number of 4.15. The 
Reynolds number range for these tests was from 1,5 x 10^ to 8.1 x 10^, 
based on free-stream conditions and the diameter of the spheres. 

Measured equilibrium-temperature distributions over the forward part 
of the sphere agreed with a laminar theory at the lower Reynolds numbers 
and with a turbulent theory at the higher Reynolds numbers for both Mach 
numbers. At a Mach number of 2.00 the recovery temperatures in the 
separated-flow region decreased slightly with increasing Reynolds number. 

Heat-transfer measurements at the stagnation point made at both Mach 
numbers agreed with laminar theory. At a Mach number of 2.00 transition 
to turbulent flow occurred at about 20° from the stagnation point. The 
heat-transfer coefficients in the turbulent boundary layer were in rea- 
sonably good agreement with a simple theory for this case. Similar 
results were obtained at a Mach number of 4.15 except that transition 
occurred farther back on the nose and, at the lower Reynolds numbers, 
the flow was laminar over the entire hemisphere. At Mach number 2.00 
the heat-transfer coefficients in the separated-flow region were about 
12 percent of the peak values on the front part of the sphere. 


INTRODUCTION 


The aerodynamic characteristics of blunt bodies at supersonic speeds 
have received considerable study in recent years. One of the principal 
reasons for this increased interest is the large reduction in local tem- 
perature that may be obtained by blunting the nose of a body. The 
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temperatures on the surface of a "blunt nose are generally less than those 
on a sharp nose at given stream conditions because of the smaller aero- 
dynamic heat-transfer rates, the increased heat-storage capacity, and the 
increased rate of heat removal by conduction in the blunt nose. 

Several experimental and theoretical investigations of the aero- 
dynamic heat-transfer characteristics of hemispherical noses or blunt 
bodies at supersonic speeds are reported in references 1 to 9* Most of 
the experimental data of these investigations were obtained at relatively 
small Reynolds numbers, and these data are generally found to be in rea- 
sonable agreement with theoretical predictions for a laminar boundary 
layer. Some measurements of the heat transfer in a turbulent boundary 
layer are reported in reference 2; however, more data are required before 
general correlations or comparisons with theoretical calculations for 
turbulent heat transfer are possible. 

The purpose of this report is to present additional experimental data 
on recovery temperatures and aerodynamic heat transfer on spherical noses 
with laminar, transitional, and turbulent boundary layers. The tests were 
made at stream Mach numbers of 2.00 and 4.15 over a Reynolds number range 
from 1.5 X 10 6 to 8.1 X 10^, based on free-stream conditions and the diam- 
eter of the spheres. The laminar and turbulent data are compared with 
appropriate theories. Heat transfer and recovery temperatures measured 
on the back of spheres in the separated-flow region are also presented. 


SYMBOLS 


C- 


’P 


pressure coefficient 


specific heat of model material 


c p 


specific heat of air at constant pressure 


D 


diameter of model 


e 


emissivity 


F 



H 


heat- transfer parameter, 
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h heat-transfer coefficient, 

k thermal conductivity 

l radius of heat-meter flange 

lh = ip + Width of air gap around plug 




T r - T w 


M 

%u 


N- 


'Pr 


P 

Q 

<1 




radius of heat-meter plug 
Mach number 

hx 

Nusselt number, - — 

K 

c v [i 

Prandtl number, — 7 — 
k 

pressure 

total heat flow rate 

heat flow rate per unit area 

Pv^D 


wall Reynolds number, 


Hw 


local wall Reynolds number, 


Pw u i x 


Pw 


stream Reynolds number. 


Pco^D 


av 


radius of model 


average radius of model wall. 


r w + r i 



temperature 

equilibrium temperature 



recovery temperature ^q w = 0 ^ 

reference temperature for radiation, taken as temperature of 
tunnel walls 

time 
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U„ 


u 


U]_* 


W 

X 

y 

5 

8 * 

e 

T) 

9 

0 * 

P 

P 

a 

T 

* 

X 


velocity ahead of how shock 
local velocity 

dimensionless velocity gradient at stagnation point. 


^1 D_ 
dx 


lx=0 


specific weight of model material 

longitudinal distance around models from stagnation point 
normal distance from surface 
boundary- layer thickness 


displacement thickness of boundary layer, 




flange thickness (see appendix) 
recovery factor, defined by equation (10) 

angular distance around models measured from stagnation point 
in spherical polar coordinate system 


momentum thickness of boundary layer, 






dynamic viscosity 
mass density 

Stefan-Boltzmann constant, 0.173 X 10"® Btu/(hr) (ft^) (°R^) 
shear stress 

axisymmetric coordinate on models in spherical polar coordinate 
system 

independent variable in heat -conduction equation for flange 


Subscripts : 
A,B 


two geometrically similar models 
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c const ant an 

f heat-meter flange 

i inside surface of model 

m material used in model 

o stagnation point 

p outside surface of heat-meter plug 

s separation point 

st stainless steel 

T location of transition 

t total 

w outside vail of model 

00 at infinity or ahead of bov shock 

1 local value just outside boundary layer (unless otherwise 

noted) 

2 reference point taken at 9 = 90 ° (unless otherwise noted) 

A bar over a symbol indicates the ratio of any temperature to stagna- 
tion temperature. 

A prime denotes values indicated by a heat meter. 

APPARATUS 


Description of Tunnels 

This investigation of pressure distributions and heat-transfer 
characteristics on spheres was conducted in three of the blowdown jets 
in the Gas Dynamics Branch of the Langley Aeronautical Laboratory. These 
jets exhaust to atmospheric pressure and are supplied with air which is 
stored in a tank field at a maximum pressure of 5,000 pounds per square 
inch and a specific humidity of less than 1 part of water per million 
parts of air by weight. The air is reduced in pressure by automatic 
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regulators to the desired stagnation pressure, which can he held con- 
stant to within l/2 percent. The stagnation temperature was also auto- 
matically regulated so that the maximum variation of this temperature 
with time was 1° F per minute. Pertinent information about the three 
blowdown jets used and the range of operating conditions available for 
the tests is given in the following table: 



M = 2 
open jet 

Variable Mach 
number closed 
jet 

M = 4 jet 

Stagnation pressure, 

lb/sq in. gage 

90 to 140 

30 to 140 

220 to 500 

Stagnation temperature, °F . . . 

120 to 680 

100 to 130 

120 to 320 

Test-section Mach number .... 

2.00 ± 0.02 

2.00 ± 0.03 

4.15 ± 0.03 

Reynolds number per foot .... 

8.5 x 10 6 

t0 6 
33.0 x 10 b 

7.0 X 10^ 
to , 

33.O x 10 b 

12.1 X 10 6 
to £ 

44.0 x 10 b 

Test-section size, 

width by height, in 

9 t>y 9 

9 t>y 9 

12 by 13 

Type of diffuser 

None 

Fixed 

Fixed 


Models and Instrumentation 

Pres sure -distribution model. - The model used to obtain the pressure 

— cr 

distributions was a - inch-diameter sphere supported by a - - inch- 

diameter sting. A sketch of this model is not shown since the shape and 
size of the model and its sting support are the same as those of the copper 
heat-transfer model shown in figure l(a) . The pressure-distribution model 
was constructed of stainless steel and was provided with 35 pressure 
orifices installed along one longitudinal line at intervals of 5° from 
0 = -10° to 9 = l60°. The diameter of the orifices was 0.02 inch. 

Mercury manometers were used for all pressures below 50 pounds per 
square inch gage, and Bourdon gages were used for pressures greater than 
this value. 

Equilibrium- temperature models . - The models used to obtain the 
equilibrium temperatures were thin-shell spheres made of Inconel. Two 
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models of different wall thickness were used. One model was 3 inches 
in diameter with a wall thickness of l/8 inch and was used to obtain 
equilibrium- temperature data at Mach number 4.15 only. This same model 

was subsequently machined to a diameter of 2-2 inches, which resulted in 

8 

a wall thickness of l/l6 inch. This latter model was then used to obtain 
additional equilibrium-temperature data at Mach number 4.15 in order to 
investigate the effect of wall thickness on the measured equilibrium 

temperatures . The 2l - inch-diameter Inconel model was also used to 

8 

measure equilibrium temperatures at Mach number 2.00. A sketch of the 
3-inch-diameter model and its sting support is shown in figure l(b) 
together with the location and method of installation of the 32 iron- 
constantan thermocouples in the sphere shell. 


The temperatures were recorded on self-balancing potentiometers. A 
manual switching system was used so that all the model temperatures as 
well as the stagnation temperature of the air could be recorded on a 
single instrument. The accuracy of the potentiometers used was l/4 per- 
cent of full-scale deflection which was 200° F for the tests at Mach 
number 2 and 330° F for the tests at Mach number 4. 


Isothermal heat-transfer model . - The isothermal heat-transfer model 
was a relatively thick-shell sphere made of electrolytic pure copper. 

The model was 3^ inches in diameter and had a 0. 3-inch- thick wall which 

remained essentially isothermal even for large heat-transfer rates because 
of the large thermal conductivity of copper. A sketch of the model, the 
sting support, and thermocouple installation is shown in figure l(a) . The 
thermocouples consisted of single constantan wires soldered at various 
locations on the model and a common heavy copper wire soldered at one 
point on the inside surface. As indicated in figure 1(a), 35 thermo- 
couples were installed at 5° intervals along one longitudinal line with 
the copper-constantan junction located approximately 0.05 inch from the 
outside surface. Also, four thermocouples were attached to the inside 
surface at intervals of 45°. The model was chromium plated to reduce 
surface abrasion during the tests. The thickness of the chromium plating 
was approximately 0.0002 inch. 

The temperature -time history of this model was obtained from a 
36-channel recording oscillograph which has elements with a sensitivity 
of 12.8 microamperes per inch of deflection. Full-scale deflection was 



An accuracy in absolute temperature of 1° F was obtained by individually 
calibrating the galvanometer elements before each test. The relative 
temperatures from any one channel were accurate to within ±l/4° F. The 
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stagnation temperature was obtained from self-balancing recording poten- 
tiometers with an accuracy of l/4 percent of full-scale deflection. Two 
instruments were used with ranges of 0° F to 600° F and 0° F to 1,200° F 
depending upon the range of the tests. 

Hemisphere-cylinder heat-transfer model .- The hemisphere-cylinder 
model was 2 inches in diameter and had a wall thickness of 0.110 inch. 

The model was made of stainless steel and was instrumented with 7 plug- 
type heat meters and 10 copper-constantan thermocouples. The location 
of the thermocouples and heat meters and the construction of the heat 
meters are shown in figure 1(c). The meters were fabricated by first 
tinning the contact surfaces of the stainless steel and constantan with 
0.0015 inch of silver solder and, then, heating the assembly under pres- 
sure to make the joints as thin as possible. The meters or plugs were 
then silver soldered into the model with a lower temperature solder and 
a final machining cut was made on the whole model. A flange that was 
0.006 inch to 0.012 inch thick held the meters in place. The air space 
was sealed at the inside surface with insulating cement and the stainless- 
steel wires were spot welded to the meter and model. 

The outside thermocouples were installed with the junction very near 
the surface. This installation was accomplished by first coating the 
bare thermocouple wires with a heat-resistant insulating paint. The 
individual wires were then inserted in separate holes drilled normal to 
the surface and spaced 0.1 inch apart, as indicated in figure l(c). The 
ends of the wires were cut off flush with the outside surface and dressed 
down smooth. The thermoelectric circuit was then completed by a very 
thin film of silver solder applied over the outside surface area including 
the two holes. 

Heat-transfer data were obtained on this model by cooling it with 
water at a constant flow rate and temperature until steady wall tempera- 
tures were indicated. The electromotive force across the constantan disk 
in the center of the heat meter is then proportional to the heat flow 
rate through the plug. Application of a correction factor to this heat 
flow rate gives the heat flow rate per unit area at the outside surface. 

The output from the heat meters was amplified by a direct-current 
amplifier which has six different scales ranging from 50 microvolts to 
2,000 microvolts full-scale deflection. The output from the amplifier 
was recorded on a self-balancing potentiometer. The overall accuracy of 
the amplifier and recorder is 1/2 percent of full-scale deflection for 
a steady input. 
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TEST PROCEDURES AND DATA REDUCTION 


Longitudinal heat conduction around the shells of the models and 
heat transfer due to radiation were included whenever necessary in the 
data reduction for the equilibrium- temperature models and the copper 
heat-transfer model. The relation used in data reduction is derived 
from the heat balance for an element of a spherical shell with small but 
finite wall thickness Ar = r w - r^ and volume r av sin 9 dSp d0 Ar where 
r av is the average or mean radius of the shell. The aerodynamic heat- 
transfer rate for an axisymmetric temperature distribution and zero heat 
transfer at the inside surface is then: 

T ref 

( 1 ) 

where in this equation T should be considered as an average temperature 
across the thickness Ar at any station 9. A similar expression for 
the aerodynamic heat transfer to a spherical shell is given in refer- 
ence 8 where no restriction is imposed on the shell thickness but the 
temperature derivatives with respect to time and 0 are assumed inde- 
pendent of r. The temperature-time derivative term in the expression 
of reference 8 is, therefore, different from the corresponding term in 
equation (l) . For the shell thickness and radius of the copper heat- 
transfer model this term in equation (l) is about 0.7 percent smaller 
than the exact value from reference 8. 

A general procedure used in all the tests was to preheat the piping 
and tunnel system up to a temperature approximating the desired stagna- 
tion temperature. This procedure insured that for the equilibrium- 
temperature tests the heat-transfer rates due to radiation were small. 

The thermal properties of air were taken from reference 10. 


q w = Wc m 


Ar - 


/ r av\ 2 8T _ k nA r /8 2 T 1 8T 


W 


St 


L w 


\d0 2 


tan 0 S9 ; 


/ 

ae V 


Equilibrium and Recovery Temperatures 

Thin-shell models .- The term "equilibrium temperature" is used to 

denote the local measured temperature when ^ = 0, and the term "recovery 

ot 

temperature" is used for the local temperature that would be obtained when 
the local value of q w = 0. Before any equilibrium-temperature data were 
obtained, the tunnel was run at steady stagnation conditions for at least 
5 minutes, and, as a result, the storage term of equation (l) was negli- 
gible. In general, when — = 0, q / 0 but has some small finite value 
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which is difficult to calculate accurately because the temperature deriv- 
atives in equation (l) are difficult to evaluate accurately from experi- 
mental data. The recovery temperatures may be obtained from the measured 
equilibrium temperatures on two thin-shell models of different wall 
thicknesses (without computing q w from eq. (l)) by plotting the equi- 
librium temperature against wall thickness and extrapolating to zero 
thickness. In general, the correct functional form of the variation of 
equilibrium temperature with wall thickness is unknown; however, a simple 
linear extrapolation to zero thickness gives the recovery temperature 
with good accuracy as is shown by the following considerations. 


If two thin-shell models, A and B, the same in every respect except 
for wall thickness and diameter, are tested in an airstream at the same 

Mach number and Reynolds number and at nearly the same temper- 


ature, then dimensional considerations require that at corresponding points 
on the two models 


''hlA _ / hlA 
V^/A VWb 


( 2 ) 


The ratio of the local heat-transfer rates would then be 


lw, A _ ^B kpo, A 1w,A ~ T r 
^w,B D A ^B T w,B " T r 


(3) 


Solving this equation for T r and substituting from equation (l) for the 
heat-transfer rates results in 


T = 
x r 


m A m 

v,B F b “ 

lA-l 

f b 


(4) 


where all temperatures are made dimensionless by dividing by the test 

F 

stagnation temperature. The ratio is given by the equation 

F B 



/ArJ 


f A 

V d O 

'a\ 

f b " 

(*sL\ 

I / 





2- 


dT 


S 2 t 


( 5 ) 


dT 


Sq 2 tan e beJ B 

^ = 0 and negligible radiation, 
ot 


for the conditions of 
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The percent error in recovery temperature due to an error in the 

... F A • 
ratio ^ is 

f B 


Tw,A _ ± 
dT r _ T w ^b 

1 


(6) 

Tr f A T Wj a 

f a 

3. W 

Ffi *w,B 

f b 



from equation (4) for given values of 

Tw,A 

and T w , b . 

Equation (6) then 

indicates that, for the conditions of 

... ..Fa 

t w,A 

** T W , B and 

f A i F B , a rather 


large error in the ratio can he tolerated without causing a corre- 

F B 


spondingly large error in recovery temperatures. Consequently, when 

D a » Db and. k^ ^ « k^ B , equation (5) may he replaced hy the approximate 

expression 


f A ^ ^A 
F b A r B 


( 7 ) 


since it follows from the conditions of T w ^a « Tw^B and the physical 
similarity of the models that the temperature derivates are also approxi- 
mately equal. Substituting equation ( 7 ) into equation (4) results in the 
approximate expression for the recovery temperature ratio 

- Ar A _ 

- ^ x w, B “ i w,A 

T-k» ^ 

^a . 1 

Ar B 


This expression shows that for thin-wall models the recovery temperature 
is given by a linear extrapolation to zero thickness of the variation in 
equilibrium temperature with wall thickness. 

The radiation term was computed for some typical tests and was found 
to be less than 4 percent of the conduction term and, hence, was generally 
neglected. The small value of the radiation term is due to the relatively 
small values of T w and T re f and, also, to the fact that these tempera- 
tures were about the same because of the preheating before each test. 
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The physical constants for the Inconel used in the data reduction 
were k m = 8.5 Btu/(hr) (ft) (°F) and e = 0.2. The models were polished 
between all tests . 


Heat-meter model .- The aerodynamic heat-transfer rates were measured 
directly on this model so that the recovery temperatures were obtained by 

plotting H (proportional to q w ) against 


T n - T, 


w 


H = 0; that is, for H = 0, q w = 0 and 


T - T 


V 


and extrapolating to 

T - T 
■*-o r 


T, 


Transient Heat-Transfer Measurements 

The aerodynamic heat transfer at Mach number 2.00 was obtained on the 
thick-wall copper sphere by a transient method. In this method the heat- 
storage or time-dependent term of equation (l) is usually the predominant 
factor in calculating the aerodynamic heat transfer. The model was pro- 
tected from the unsteady starting conditions in the jet by disposable 
water-cooled covers. The covers consisted of double-walled hemispherical 
shells which were fitted together with a leak-proof seal. After attaining 
the desired steady conditions in the jet, a quick-release mechanism was 
actuated and the covers were blown off the model by the airstream. In 
this way the model was maintained at a uniform temperature until the 
covers were released. The data indicated that this initial temperature 
distribution was uniform to within 2° F over the entire region where 
temperatures were measured. 

The temperature variation with time at two points on the copper sphere 
during a typical test is shown in figure 2(a) . The complete temperature- 
time history up to 1 second from time zero for the same test is shown in 
figure 2(b). The initial temperature of the model for this test was 85 ° F 
as indicated in the figures. The outside temperatures were measured 
approximately 0.05 inch from the outside surface (fig. 1 (a)) and the 
inside temperatures, denoted by the flagged symbols in figure 2 (b), were 
measured at the inside surface. 

Figure 2(a) shows that the temperature derivatives with respect to 
time at small values of t depend on the radial location. However, in 
the derivation of equation (l), it is assumed implicitly that the temper- 
ature derivatives are independent of the radial location. Consequently, 
any temperature data for 0 <t< 0.3 were not used. 

Figures 2(a) and 2(b) show that the temperature also varies consider- 
ably with radial location. This variation indicates that small errors in 
the depthwise location of the outside thermocouples would cause errors in 
the derivatives of T with respect to 9 required in equation (l). These 
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errors were minimized by adjusting the time zero for each thermocouple 
record according to the radial location of the thermocouple. This adjust- 
ment was obtained by disregarding the early part of the records (gener- 
ally, t < 0.5 second) and extending the temperature-time curves backwards 
to the initial temperature by using the same slope and curvature that 
was found at the larger times (generally, t > 0.3 second). This pro- 
cedure is illustrated in figure 2(a) where the extended portions of the 
curves are shown as dashed lines. The intersection of each extended curve 
with the line for the initial temperature of the model was then taken as 
the adjusted n time zero" for that particular thermocouple location (this 
adjusted time zero in fig. 2(a) is at -0.12 second for the outside thermo- 
couple and 0.l6 second for the inside thermocouple), and the temperatures 
were then read again at the adjusted times. The result of using this pro- 
cedure on the typical test of figures 2(a) and 2(b) is shown in figure 2(c). 
The adjusted temperatures derived from the inside thermocouples are now 
about the same as the adjusted outside temperatures and some of the appar- 
ent irregularities in the temperature distributions of figure 2(b) have 
been removed. This agreement indicates that the general procedure just 
described is valid, and plots of the type of figure 2(c) were then used 
to evaluate the derivatives of T with respect to 0 required in equa- 
tion (l). No data at t > 0.6 second were used from this model because 
of the increasing ratio of the conduction to storage terms at the larger 
times. This limitation minimized the errors involved in evaluating the 
conduction since the ratio of the conduction to storage terms was always 
less than about 0.3 for t < 0.6 second. The heat transfer due to radia- 
tion was found to be negligible as compared with the storage and conduc- 
tion terms. 

The physical constants for copper used in the data reduction were 
W = 559 Mj/cu ft, c m = 0.0915 Btu/(lb)(°F), k m = 220 Btu/(hr) (ft) (°F) , 
and e = 0.5. The model was polished before each test. 


Direct Measurement of Heat Transfer 

A plug- type heat meter was used for the direct measurement of aero- 
dynamic heat-transfer rates on the hemisphere-cylinder model which was 
tested at Mach number 4.15. A general objection to the plug- type heat 
meter has been that the plug may cause large perturbations in the local 
wall temperatures of the model. These temperature perturbations in turn 
cause disturbances in the local boundary-layer characteristics and, also, 
may introduce extraneous heat-conduction effects in the plug itself. The 
heat meter and test procedure used in this investigation were designed to 
minimize these effects. The model was cooled with large amounts of water 
(about 100 pounds per minute) in the manner indicated in figure 1(c) and, 
thereby, a nearly uniform temperature was maintained at the inside sur- 
face of the plugs and surrounding model wall. The outside temperatures 
of the plug and surrounding model were also about the same because of the 
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relatively thin wall and because the materials used in the plug (constan- 
tan and stainless steel) were chosen so that the net thermal conductivity 
through the plug and model wall was nearly the same. The overall effect 
of the thin wall and large cooling rates was to cause most of the heat 
flow to take place in a direction normal to the wall in the model as well 
as the plug. The air space around the plug (see fig. 1(c)) forces a one- 
dimensional heat flow through the constantan disk. After steady temper- 
atures have been attained, the total heat flow rate per unit area through 
the plug is then 


<4 = 



(8) 


where £T C is the temperature difference across the constantan disk and 
Ar c is the thickness of the constantan disk. This temperature difference 
is obtained from the electromotive force appearing across the stainless- 
steel wires to the heat meter. (See fig. 1(c).) The thermoelectric out- 
put of the stainless steel and constantan junctions is known as a function 
of temperature from a previous calibration of the same materials. 

The aerodynamic heat-transfer coefficient is given by the equation 


h - /h\ 

T r - T p ' \h '/ 


( 9 ) 


where — is a factor that corrects for the amount of heat transferred 
h' 

between the surface of the flange and the airstream and, also, for the 
amount of heat conducted into the flange from the model. Heat transfer 
due to radiation is not included in equation (9) since it was negligible 
in comparison with the aerodynamic heat transfer. The surface tempera- 
ture of the plug Tp' is computed by assuming one -dimensional heat flow 
through the plug and by using the temperature measured at the inside 
surface. The temperature distribution along the surface of the flange 

determines the correction factor -jik and depends on the heat-transfer 

h 

coefficient and the dimensions and thermal conductivity of the flange as 
well as the surface temperatures of the plug and model. The temperature 
distribution on the flange may be calculated from a differential equation 
which is derived in the appendix. 

The data presented in this report have not been corrected for the 
flange effect because of possible experimental errors in the measured 
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temperatures T w and T p ; however, typical values of the correction 
factor and the effect of various parameters on the factor are given in 
the appendix. 

Values of T p , T w , and are shown for a typical test in fig- 
ure 2(d) . The temperatures ^w and axe measured on the outside 

and inside surfaces of the model at the locations indicated in figure l(c). 
Possible errors in these measured temperatures and the derivation of 
equation ( 9 ) are also given in the appendix. 

The physical constants for the constantan and stainless steel used for 
this model were k c = 12.8 Btu/(hr) (ft) (°F) and k s t = 9*5 Btu/(hr) (ft) (°F) , 
respectively. The stainless steel was type 3 O 3 , that is, about 18 per- 
cent chromium and 9 percent nickel. The thermoelectric output of the 
particular stainless-steel — constantan combination used in this model is 
given in the following table: 


Average temperature of 
junctions, °F 

Microvolts per 

Op 

60 

27.5 

80 

25.4 

100 

24.2 

120 

25.4 

140 

23.1 


This model was polished before most of the tests; however, several tests 
were made without polishing in order to observe the effect on the heat 
transfer of the natural- surface abrasion caused by foreign material in 
the airstream. 


RESULTS AND DISCUSSION 


Pressure Distributions 


The variation of pressure coefficient Cp with 0 at M = 2.00 for 
various stream Reynolds numbers is shown in figure 3» Over the front part 
of the sphere the pressure does not vary appreciably with Reynolds number 
and up to about 0 = 60° is closely approximated by the Newtonian 
distribution 


c p = 


■'p,o 


cos^9 


1 6 


MCA TN 4125 


where C p 0 was calculated from the average Mach number of the jet and 
the normal shock relations. On the hack of the sphere within the 
separated-flow region the Reynolds number has a large effect on the pres- 
sure and also on the separation point. 

A comparison of the separation point obtained from schlieren photo- 
graphs of the recovery-temperature model with the separation point indi- 
cated by the pressure data of figure 5 is shown in figure 4. In general, 
the flow separates somewhat farther forward on the recovery- temperature 
model, as shown by the schlieren data, than on the pres sure -distribution 
model. The different location of the separation points on the two models 
is probably caused by the different sting diameters or some other change 
in the flow associated with the particular jets used for the tests. The 
recovery-temperature model has a ratio of sting diameter to model diameter 
of 0.24 as compared with a ratio of 0 .l 8 for the pressure-distribution 
model. (See figs. l(b) and 1(a).) The schlieren photographs were taken 
of the recovery-temperature model in the variable Mach number closed jet 
and the pressure data were obtained in the M = 2 open jet with the nose 
of the model 7/8 inch inside the end of the jet. Sample schlieren photo- 
graphs at three Reynolds numbers are shown as figure 5 to illustrate the 
rearward movement of the separation point with increasing Reynolds number. 
The location of separation was assumed to be at the point of intersection 
of the forward oblique shock with the surface of the sphere. 

The pressure distributions obtained on the sphere at M^ = 4.15 and 

at Roo = 6.7 X 10^ and 9.1 X 10^ are shown in figure 6 . The pressure 
over the forward part of the model is in good agreement with that calcu- 
lated by the Newtonian expression almost all the way to the separation 
point. This small change in Reynolds number had little effect on the 
separation pressure or location. 


Equilibrium and Recovery Temperatures 

^ = 4.15.- The thin-wall Inconel model was tested first with a 
wall thickness of l /8 inch at = 4.15. After these tests were com- 
pleted, the same model was machined to a wall thickness of l/l 6 inch and 
more data were obtained. The results of both series of tests are pre- 
sented in figure 7 where the ratios of the measured equilibrium tempera- 
tures to the stagnation temperatures T e /T 0 are P lQ tted against stream 
Reynolds number R^ for values of 9 from -10° to 145°. On the front 
part of the sphere from the stagnation point to about 9 = 35 ° > there is 
little dependence of T e /T 0 on R^ (fig. 7(a)). From about 9 = 35° 
to 9 = 115° (figs. 7(b) to 7(d)), T e /T 0 generally increases with 
Reynolds number R^. For 9 > 11 5° (figs. 7(a) and 7 (e)), T e /T 0 is 

about constant when R^ > 4.5 x 10^. The values of T e /T 0 are higher 
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on the thin-wall model than on the thick-wall model at all values of 9 
except in the range of 9 from approximately 90° 1° 115° • This tempera- 
ture reversal is caused by the increased heat conduction into this region 
on the thick-wall model since the minimum temperatures are found here. 

Transition occurred at Roo 158 4.5 X 10^ at all values of 9 > 45° 
as indicated by the abrupt increase in T e /T 0 on the thin-wall model at 
this Reynolds number. For values of 9 > 115° and R^, < 4.0 x 10 , the 
measured temperatures on the thick-wall model may be about 1 percent too 
low because of insufficient testing time to allow for the extremely low 
heating rates in this region. 

The data shown in figure 7 have been extrapolated linearly to zero 
wall thickness since, according to the discussion in the section entitled 
"Thin-shell models," this procedure is permissible if T w ^ 83 T w,B and 
D A » Dp. The results of the extrapolation are shown in figure 8 where 
the ratio of the recovery temperature to the stagnation temperature T r /T 0 
is plotted against 9 for three Reynolds numbers. The curves labeled 
"laminar theory" and "turbulent theory" were computed from the relation 



( 10 ) 


where r| = ^/Np r for the laminar-flow theory, t} = ^/Np r for the 

turbulent -flow theory, and the Prandtl number was assumed constant at 
Np r = 0.7. Just as is shown in figure 7 , these data are almost independ- 
ent of R^ for 9 < 40°, and at = 3-0 X 10° the data are slightly 
below the laminar -theory curve. For 9 > 40°, T r /T 0 at the larger 
Reynolds numbers is higher than the laminar data and tends to follow the 
trend of the turbulent-theory curve. This behavior indicates that transi- 
tion occurred at about 9 = 40° to 45°. Downstream of the separation 
point which was between 9 = 90° and 105° (in agreement with the pressure 
data of fig. ( 6 )) T r /T 0 is again independent of R^ and increases with 
increasing distance around the model. All the data up to 9 = 9 O 0 would 
be in better agreement with the theory if the laminar -theory and turbulent- 
theory curves were reduced about 1 percent in order to agree with the data 
at 9 = 0°. 

Some recovery-temperature data from the heat-meter model are also 
shown in figure 8 . These data were obtained from figure 9 where the 
dimensionless heat-transfer parameter H is plotted against the tempera- 
T - T 

ture parameter -2- — . In accordance with a procedure given in the 

T o 
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section entitled "Heat-meter model, " the recovery temperatures from this 
figure are computed from the intersection of the faired curves with the 
axis H = 0 so that the expression for T r /T 0 is 




H=0 


Note, however, that the surface temperature distribution on the heat-meter 
model is such that the condition q w = 0 is generally found at only one 
point or area on the model during any one test. Thus, in figure 9 the 

( T - T \ 
smallest values of —2^ — ) were all obtained 

in the same test, so that for this test q w « 0 at 9 = 75° and- 90°; 
for 9 < 75° the wall temperatures were still considerably below recovery 
temperatures. According to a qualitative analysis of the effect of wall 
temperature distribution on the local recovery temperatures, the recovery 
temperatures from figure 9 should then be somewhat less than the values 
on the thin- wall spheres where for any one test q w « 0 over the entire 

model. Comparison of the data in figure 8, however, shows good agreement 
between the two sets of data at comparable Reynolds numbers. Apparently, 
then, this effect is too small to be measured in the present tests. 

= 2.00.- The ratio of the equilibrium temperature to the stagna- 
tion temperature T e /T 0 measured at a stream Mach number of 2.00 is pre- 
sented in figure 10. These data were all obtained on the Inconel model 
with a wall thickness of l/l6 inch. On the front part of the sphere the 
results are similar to the recovery-temperature data at = 4.15 

(fig. 8) except that for the larger Reynolds numbers ( R^, ^ 6.45 X 10^) 
transition apparently occurred somewhat farther forward at 9 « 25° to 50°. 

At the lower Reynolds numbers ( R^ S 5*49 X 10^) transition apparently 
occurs ahead of 9 = 90° as indicated by the increasing values of T e /T 0 
in the vicinity of 9 = 70° to 80°. 

The theoretical curves were computed by use of equation (10) with 
the same values of T). The agreement between the data and the theory 
would be improved if the theoretical values were reduced by about 
l/2 percent. 

In the separated region (fig. 10) , T e /T Q was essentially constant 
with increasing distance around the model but decreased slightly with 
increasing Reynolds number. This decrease is probably caused by the same 
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flow mechanism which caused the corresponding decrease in the pressure 
in the separated region shown in figure 5» 


Heat-Transfer Coefficients 

Stagnation point . - All the heat-transfer data obtained at the stag- 
nation point in this investigation are shown in figure 11 as a "wall" 

Nusselt number — plotted against a parameter which is the product of 
k w pTJ D 

a "wall" Reynolds number Rp = — - and the dimensionless velocity 

gradient up* evaluated at the stagnation point. The recovery tempera- 
tures used to compute h were taken from the data of figures 8 and 10. 

The heat-transfer data at = 2.0 were obtained on the thick-wall 
copper sphere by the transient method, and the data at ^ = 4.15 were 
obtained on the hemisphere-cylinder model with the heat meters. The heat- 
meter data are shown uncorrected for any errors caused by the heat trans- 
ferred to the flange. Some of the data at = 4.15 were obtained with 
the hemisphere-cylinder model yawed 15 ° in order to provide a cross check 
of the heat-transfer coefficients obtained from the different meters. 
Comparison of these data shows that the heat meter located at 9 = 15° 
in the unyawed position and at 9=0° when the model was yawed (see 
fig. 1(c)) was indicating heat-transfer rates about 10 percent higher than 
the meter at 9=0° (model unyawed) . This discrepancy is probably caused 
by constructional variations in the plugs since yawing the model should 
have no effect on the heat transfer on the hemispherical nose. 

Lines shown in figure 11 labeled "laminar theory" are computed from 
the results of reference 9 where it is shown that in the vicinity of the 
stagnation point the local Nusselt number divided by the square root of 
the local Reynolds number or the quantity 


\ 


is a function only of the Prandtl number and the ratio of vail tempera- 
ture to stagnation temperature. If the dimensionless velocity gradient 
at the stagnation point u-^* is introduced into this quantity, then 
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The values used for — — were O.663 and 0.639> respectively, for 

T w /t q = 1.0 and 0.5 (ref. 9). The values used for u^* were 1.54 for 
l^o = 2.00 and 1.19 for 1^, = 4.15, as obtained from the pressure- 
distribution data of figures 3 and 6. 

The data at = 2.00 are in good agreement with the theory but 
the data at = 4.15 are about 12 percent higher than the theory. 
Application of a correction for flange effects would reduce these data 
at ^ = 4.15 by 15 to 25 percent; this reduction depends on the plug 
and wall temperatures and flange thickness, as discussed in the appendix. 
Apparently, then, the corrected data at ^ = 4.15 would be somewhat 
lower than the theory by an amount which is well within the uncertainties 
due to flange correction and constructional variations in the plugs. No 
effect of a variation in T w /T 0 as indicated by the theory could be 
detected in the data. 

Variation of heat-transfer coefficients with 9 at = 2.00.- 
The heat-transfer coefficients on the copper sphere at M^, = 2.00 are 
presented in figure 12 as the ratio of the local heat-transfer coefficient 
to the value at the stagnation point. The local heat-transfer rates were 
calculated from the experimental temperature -time histories by using equa- 
tion (l) , and the recovery temperatures were obtained from the experi- 
mental data of figure 10. 

The large increases in heat transfer occurring on the front part of 
the sphere and reaching a maximum at about 9 = 40° (fig. 12) are evi- 
dently caused by transition to turbulent flow, since the heat-transfer 
coefficients at the stagnation point are in agreement with laminar theory 
(fig. 11). Also, if the boundary layer were entirely laminar, the heat 
transfer would be expected to decrease with increasing 9 as shown in 
previous investigations. (See, for example, refs. 4 and 8.) 

The heat-transfer distribution for a laminar boundary layer, as shown 
in figure 12, was computed from the method of reference 11 by use of 
Mangier 's transformation (ref. 12). Comparison of the theory with the 
experimental data shows reasonably good agreement from the stagnation 
point up to the region of 9 = 10° to 30° where large increases in h 
occurred. This value of 9 is presumably the location for the beginning 
of transition and is subsequently referred to as 9ip. 
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The theoretical heat-transfer distribution in a turbulent boundary 
layer was computed by using a modification of Falkner's expression 
(ref. 15 ) for the skin friction on a flat plate given by the equation 


'w 


and Reynolds analogy in the form 


0.0131 


PVu/ /(yuix U/7 


Hv / 


Pw u l c p 


= 1.2 


'w 


Pw U l 


where the gas properties are evaluated at the local wall temperature and 
local pressure. (Skin-friction coefficients and Reynolds analogy from 
these formulas are compared with recent experimental data on flat plates 
at supersonic speeds in refs. 14 and 15 . ) The resulting expression for 
the heat-transfer coefficient is 


h = 0.0157 


Pw u l c p 



( 12 ) 


Combining equation (12) with the expression for the laminar heat-transfer 
coefficient at the stagnation point obtained from equation (ll) and 
assuming constant T w gives the equation 


h 

h e=o 


0.0157 Np r 



' p i u i \ 6 / 7 



'p U D 

W °o 


5/14 


x=0 


(13) 


which was used to compute the long-short-dash curves in figure 12 labeled 
’’local flat plate (turbulent theory).” The local values of pressure and 
velocity were obtained from the measured pressure distributions of fig- 
ure 3 and x was measured from the stagnation point. The two curves 

shown were computed for Rj) = 11 x 10 ° and 19 X 10 which are the 
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approximate limits of this parameter for the tests. Comparison of the 
experimental data with the theory indicates that equation (lj) predicts 
the correct variation of the heat transfer with 9 from 9 = 40° to 90° • 
The predicted magnitude of the heat transfer on this portion of the 
model is generally conservative. The experimental values are as much 
as 20 percent below the computed values except for the test at 

= 2.74 X 10^ where the data are as much as 40 percent below the 
theory. The agreement between the experimental data and the theory 
indicates that, within the experimental scatter of the present data, 
the simple flat-plate formulas given by equation (12) predict the tur- 
bulent heat transfer on spheres with reasonable accuracy. No consistent 
trend caused by changes in Rp or T w /T 0 is evident in the data. The 
experimental scatter in the data is probably masking any such effects 
over the small range of Reynolds number and temperature available in 
these tests. The heat-transfer coefficients in the separated region 
on the back of the sphere are about 12 percent of the peak turbulent 
values . 

Heat-transfer coefficients at = 4.15*- The experimental heat- 
transfer coefficients on the hemisphere-cylinder model at = 4.15 
are presented in figure 13 in the form of the wall Nusselt number divided 
by the square root of the wall Reynolds number. This particular param- 
eter was used here rather than h/hg_Q in order to correlate the laminar 
heat-transfer data and, also, to facilitate comparison with other experi- 
mental and theoretical investigations, such as those of references 4 
and 9* The heat-transfer rates were obtained directly from the heat- 
meter readings by means of equation (8) with no correction included for 
flange effects. (See the appendix for a discussion of these corrections.) 
The recovery temperatures were obtained from the data of figure 8 and 
the local flow quantities were computed from the pressure distributions 
of figure 6 and the measured wall temperatures . Each data point shown 
in figure 13 is the arithmetic average of the values obtained from the 
number of tests indicated in the key. 

The curves shown in figure 13 for the heat transfer in a laminar 
boundary layer were obtained by computing the variation in heat transfer 
around the sphere by the method of reference 11 and applying a correction 
to these results in order to bring the stagnation-point value into agree- 
ment with the theory of reference 9. The local flow quantities needed 
in this calculation were computed from the experimental pressure- 
distribution data of figure 6. 
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The theore tical distribution of the laminar heat-transfer parameter 
hx/k w /pwU]_x/p w in a turbulent boundary layer is obtained from equation (12) 
and is given by the relation 




hx 

P^l x 
l J -w 


= (0.0157) Npr 


15/14 5 / 1 ^ 

u lPl x / p v u i x \ 

(uiPlx) 2 \ hr ! 2 


for T w = T w ^ 2* Th e subscript 2 denotes quantities evaluated at any 
convenient reference point which in this case is taken at 0 = 90° • 

The data of figure 13 have been divided into two sets according to 
the condition of the model surface during the tests. One set of data was 
obtained when the model was carefully polished before each test. These 
data are designated by the open symbols in figure 13 and are considered 
representative of a "smooth" surface. The other set of data, designated 
by the solid symbols, was obtained when no attempt was made to keep the 
model polished and are considered representative of a "rough" surface. 

This latter procedure resulted in a considerably rougher model surface 
than the former procedure because of the cumulative effect during the 
tests of the natural abrasion caused by foreign material in the airstream. 

In figure 13 a comparison of the data for a smooth surface (model 
polished before each test) with the laminar theory indicates that, when 

R^ = 3.3 X 10^, the data are about 15 to 30 percent higher than the theory 
but tend to decrease with increasing 0 by about the same amount as the 

theory. For = 4.9 x 10^, however, the values of the heat-transfer 
parameter at 0 = 45°, 60°, and 75° are from 50 to 200 percent greater 
than the laminar theory and tend to approach the level predicted by the 
turbulent theory. Apparently, transition from laminar to turbulent 
boundary layer occurred between 0 = 30° and 45° at the larger Reynolds 
number. Note that alternate meters are located on opposite sides of the 
model (fig. 1(c)) so that if turbulent flow existed on one side of the 
model only, then alternate meters would indicate higher heat transfer. 

The data for the rough surface (model unpolished) are, in general, 
higher than the other data in figure 13 except at 0=0° and 15° • At 
0 = 30°, 45°, and 75 0 for the lowest Reynolds number, these data are 
about 15 percent higher than the corresponding data from the smooth sur- 
face. As the Reynolds number is increased, the heat transfer tends to 
get progressively larger for all values of 0 = 30° and approaches the 
levels predicted by the turbulent theory. The data at 0 = 90° are, in 
general, below the turbulent theory. This may be attributed to two factors. 
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One of the factors is that the meter at 9 — 90° no ^ exposed to ero- 
sion effects to the same extent as the other meters and, hence, the local 
surface would not be as rough. Also, the small measured values of the 
heat-transfer rates at 9 = ^^ } as shown for typical conditions in fig- 
ure 9, result in larger errors in the heat-transfer coefficient at this 
location since the instrumentation errors, as discussed in the appendix, 
would have a larger effect. 

Application of a flange correction (see the appendix) would reduce 
all these data by 15 to 25 percent. This correction would bring the 
smooth-surface data at the smaller Reynolds numbers into good agreement 
with the laminar theory. The higher values of the heat transfer on the 
rough surface would be within the upper and lower limits of the turbulent 
theory. Evidently all other values of the heat transfer are in the 
transitional range. If the larger experimental values are for a fully 
turbulent boundary layer, then the flat-plate formulas again predict con- 
servative values of the heat transfer. 

Typical values of heat-transfer rates obtained on the hemisphere- 
cylinder model at R^, = 3*3 x 10^ and H» = 4.15 are shown in figure 9* 
In order to correlate the data for a laminar boundary layer, the dimen- 
sionless parameter H is plotted against the temperature parameter 

~ The slope of the faired lines would be the local stream Nusselt 

T 0 

number divided by the square root of the stream Reynolds number, and the 

T 0 - T r 

value of the temperature parameter when q w = 0 would be . The 

relative scatter at 9 = 60°, 75°, and 90° is greater than at 9=0°, 
15°> 30°, and 45°. This is probably caused by the greater sensitivity 
of the flow to surface roughness at the rearward stations. The heat- 
transfer-coefficient data presented in figures 11 and 13 were not derived 
from plots of the type shown in figure 9 since a larger number of tests 
at a given Reynolds number is then required in order to evaluate the 
Nusselt number. 


Effect of Roughness on Transition 

The quantitative effect of roughness on transition cannot be deter- 
mined from the present tests since no attempt was made to measure or con- 
trol the actual roughness. Some qualitative information may be gained 
from the tests, however, by calculating displacement and momentum thick- 
nesses for comparison with the possible roughness present. Also, momentum- 
thickness Reynolds numbers at transition may be useful for comparison with 
other data. 
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The displacement and momentum thicknesses of the laminar boundary 
layer on a sphere have been calculated by the method of reference l6 by 
use of the experimental pressure distributions of figures 5 and 6. The 
displacement thickness times the square root of the wall Reynolds number 
was found to be almost constant over the first 50° of the sphere; that 
is, 


6 * 

D 





0.5 


for 0 < 30° and for the range of Mach numbers and wall temperatures of 
the present tests. For the maximum Reynolds number range and the diam- 
eters applying to the heat-transfer tests at Mach numbers 2.00 and 4.15, 
the minimum value of 5* is about 0.0003 inch and the maximum value is 
about 0.0005 inch on both sets of tests. These small values indicate 
that very small roughness heights of the order of 0.0002 inch would reduce 
considerably the transition Reynolds numbers if correlations obtained at 
low speed (ref. 17) of the effect of roughness can be applied to the 
subsonic flow region on the sphere. Observations made of the abrasion 
damage to the surface of the various models indicated that, in spite of 
polishing before the tests, at least this amount of roughness was present 
after the tests. The relative damage to the models was, in general, the 
greatest on the copper model and the least on the Inconel recovery- 
temperature models. 

The Reynolds numbers for transition based on the local momentum 
thickness 0* and the local flow quantities at the edge of the boundary 
layer are plotted against angular distance around the model in figure 14 
for all the data obtained in this investigation. The location of transi- 
tion was arbitrarily taken as the point where the heat-transfer coefficient 
or recovery temperature first began to increase from the nominal laminar 
values. Figure 14 shows that most of the data from the heat-transfer 
models with the rough surface (model unpolished) at M^, = 2.00 and at 
M^ = 4.15 are lower than the rest of the data and correlate approximately 
on one line. For these particular tests, this line apparently represents 
a minimum critical Reynolds number below which the flow was laminar. The 
data above this line are, in general, from the polished-model data at 
^ = 4.15 and from the data for the recovery-temperature models. The 
two points plotted at 0 = 90° are from the smooth- surface model for 
Roo ^ 3.3 X 10^ and are included to show that transition occurred down- 
stream of 0 = 90° for these conditions. All the data indicate that 
Rq* T tends to increase with increasing 0. 
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CONCLUDING REMARKS 


An investigation was made of the pressure and equilibrium-temperature 
distributions on a sphere at Mach numbers of 2.00 and 4.15- The local 
aerodynamic heat transfer was also measured on a sphere at Mach num- 
ber 2.00 and on a hemisphere-cylinder at Mach number 4.15. The maximum 
Reynolds number range for these tests was from 1.5 k 10^ to 8.1 x 10 , 
based on free-stream conditions and the diameter of the sphere. 


The pressure distributions on the front part of the sphere were 
independent of the Reynolds number and in good agreement with the theo- 
retical Newtonian distributions at both Mach numbers. At the lower Mach 
n um ber the location of separation and the pressures in the separated-flow 
region were affected by the Reynolds number. 

The equilibrium temperatures on the front part of the sphere at both 
Mach numbers were about 1 percent lower than those predicted by a simple 
theory which assumes that the local recovery factor (defined in terms of 
the local static temperature) is the square root of the Prandtl number 
for laminar flow and the cube root of the Prandtl number for turbulent 
flow. The data at the lower Reynolds numbers agreed with the trend given 
by the theory for laminar flow. As the Reynolds number was increased, 
the temperatures tended to increase and agreed more closely with the 
theory for turbulent flow. 

Downstream from the separation point the recovery temperatures for 
the tests at Mach number 4.15 increased with increasing distance around 
the sphere and were practically independent of Reynolds number . At Mach 
number 2.00 the equilibrium temperatures were more nearly constant in the 
separated-flow region and decreased slightly with increasing Reynolds 
number . 

The heat-transfer coefficients at the stagnation point from the 
tests at both Mach numbers were in good agreement with the theory of 
Reshotko and Cohen (NACA Technical Note 3513)- The velocity gradients 
used in the theory were evaluated from the experimental pressure data. 

The local heat-transfer coefficients from the tests at Mach num- 
ber 2.00 were in agreement with the theory for laminar flow in the region 
up to approximately 20° from the stagnation point. Large increases in 
heat-transfer coefficient caused by transition to turbulent flow occurred 
beyond this region with the maximum values observed at about 40° from the 
stagnation point. From this location back to the separation point the 
approximate level in the heat-transfer coefficient as well as the varia- 
tion with distance was fairly well predicted by a simple theory based on 
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turbulent heat-transfer formulas for flat plates. The heat-transfer 
coefficients in the separated-flow region were about 12 percent of the 
peak turbulent values on the front part of the sphere. 

The heat-transfer coefficients from the tests at Mach number 4.15 
at the higher Reynolds numbers were generally similar to the results at 
Mach number 2.00, except that transition usually occurred farther back 
on the nose. The maximum values of the heat -transfer coefficient down- 
stream of transition were again in reasonably good agreement with the 
simple turbulent -flow theory. For Reynolds numbers of about 5*0 x 10 
or less, the data over the entire hemispherical nose were in agreement 
with the theory for laminar flow. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., August 6, 1957* 
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APPENDIX 

HEAT-METER CORRECTIONS 


The largest error in the heat-meter data is caused by heat conduction 
in the metal flange which bridges the gap between the plug and the model 
(fig. 1(c)). The external surface of this flange is exposed to approxi- 
mately the same unit heat-transfer rates as the plug and model surface. 

The insulating properties of the air gap force this heat to flow into the 
plug or the surrounding model. The actual quantity of heat conducted 
into the plug from the flange depends on the temperature distribution in 
the flange and the temperatures of the plug and model. The indicated heat- 
transfer rate as obtained directly from the temperature drop across the 
constantan disk is then corrected by an amount which depends essentially 
on the difference in the actual surface area involved in the total heat 
flux through the plug and the projected area of the disk. 


Ratio of Corrected to Indicated Heat-Transfer Coefficients 

A section of a heat meter with symbolic notations is shown in 
sketch 1: 



Sketch 1 

The total heat fy. flowing per unit time through the constantan disk of 
thickness Ar c is the sum of Q,p into the surface area of the plug 
(the projected area of the constantan disk) and conducted in from 

the flange; that is, 


Qt - Qp + Of 


(Al) 


NACA TN 4125 


29 


where 


Qp - h(T r - Tp^itZp ► 

^ = k st(|l) 2 *V 

t-tp 


(A2) 


The temperature gradient — — is assumed constant across the flange 
thickness e. 


An "indicated" heat-transfer coefficient is defined as 

AT„ 


h' = 


Ar, 


T - T ' 
1 r x p 


(A5) 


where T 1 is the surface temperature of the plug, which is calculated by 
assuming P one-dimensional heat flow through the entire plug. The ratio of 
the indicated heat-transfer coefficient to the true value is then obtained 
by combining equations (Al) , (A2) , and (A3) . This ratio is 


h4 

h 




_d_ 

dZ 



l=l v 


(A4) 


where 


*P 


st 


. 1 P 


From the notation shown in sketch 1 and 
dimensional heat flow, the value of T' is 


the assumption of one- 
given by the relation 


P 


+ 


k c Ar X + Ar 2 \ 

^ A-c / AC 


+ Ti 


,(A5) 
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where £T C = T 2 - Tq. A more accurate value of Tp may be computed from 
equation (A4) by assuming that 


%> ^ k st( T p ^ 2 ) 

(T r - T p )„I p 2 ar 2 (T r - T p ) 


(a6) 


and by substituting equation (A 5 ) for h'. 


The result is 


T n = T 2 + — — “ 

P 2 k st ^c C * 


(A7) 


where 


A 




_d_ 

dl 



1=1 


P 


and from sketch 1 


t 2 = 



+ T i 


(A8) 


(A9) 


Subtracting equation (A 7 ) from equation (A5) and using equation (A 9 ) then 
gives 


V 



(A10) 


Adding ^T r - Tp') to both sides of equation (A10) and using equation (A 3 ) 
for £T C results in the expression 


T, 




(All) 


NACA TN 4125 


31 


The correction factor is calculated from equations (A4) and (All) 

h 


for known values of h' and A. An_iterative procedure is necessary 
since A is a function of h and 




Temperature Distribution on the Flange 
dT 


-==- at the inside edge of the flange 
di 


The temperature gradient 

depends on the quantity of heat transferred to the flange from the air- 
stream and the surrounding model as well as the thermal conductivity of 
the flange and its dimensions. These effects may be approximately 
accounted for by writing the heat balance for an element of the flange 
of width dZ and assuming a constant temperature across the thickness e 
of the flange. The heat-transfer coefficient h and the recovery temper- 
ature T r are assumed constant throughout. The resulting differential 
equation is 


(T_ - 
r 


+ I A.h - t') - 

Z dll r ) 


st 


e ( T r 


- T) = 0 


(A12) 


with the boundary conditions of 


Z = Z p 

z = z h 


T = T t 


T = T. 


w 


(A13) 


Introducing the independent variable 


X = 


k st € 


into equation (A12) results in 


- T ) 4 Sc( T ^ - T) - (T r - T) = 0 


(A14) 


dX‘ 
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which has the general solution 

= 1 - Cii 0 (x) - C 2 K 0 (X) 
1 r 


(A15) 


where Iq and Kq are modified Bessel functions of the first and second 
kind of order zero and Cp and C 2 are arbitrary constants chosen to 
satisfy the boundary conditions (eqs. (AI 3 )). Note that these boundary 
conditions imply a discontinuity in temperature gradient at the edges of 

riT 

the flange since — = 0 for Z_ > l > l h and no boundary conditions 
dZ P n 

for the values of — on the flange are imposed. It is believed, there- 

dZ 

fore, that the results obtained by using the boundary conditions 
(eqs. ( A 13 ) ) are conservative since the physical requirement of a con- 
tinuous temperature gradient would tend to reduce the temperature gradient 
at the inside edge of the flange and, thereby, to decrease the quantity 
of heat conducted into the plug from the flange. 

The temperature gradient as obtained from equation (AI 5 ) is 


_d_ 

dZ 



Z=Z p 



- c i i i( x p) 


* C 2 K l(Xp) 


(A16) 


where Ip and Kp are modified Bessel functions of the first and second 
kind of order one. The quantity A then becomes 


A = 1 + 


*P 


li.h' 


-C-,1 


1 X 1 


(*p) + C 2 K l( X p) 


(AIT) 


The reciprocal of A is plotted against 
values of Xp. 


T v - T- 


Values of h', Xp' = 


T r - T p 


^w - 


£ in figure 15 for typical 


^ ' Z„, and — — from the test shown 

ks t e P’ T r - T p ' 


in figure 2(d) are given in the following table: 
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9, 

deg 

© 

in. 

© 

h 1 , 

Btu 

© 

S' 

© 

** - V 

© 

1 

X 

(Fig. 15 
and 

col. (5)) 

© 

Tr - T p 

© 

h 

h' 

(Eq. 

(A18)) 

© 

v - v 

© 

h 
h ’ 

T - T 1 

r / p 
(Eq. 

(All)) 

T r - Tp' 

T r - T p ' 

(sec) («2) (Op) 

15 

0.006 

0.0543 

1.068 

0.034 

O.63 

1.018 

0.619 

0.018 

0.71 

15 

.003 

.0543 

1.510 

.034 

• 75 

1.013 

.720 

.018 

.78 

45 

.006 

• 0375 

.887 

.024 

.63 

1.012 

.623 

.013 

.71 

75 

.006 

.0144 

.549 

.015 

• 55 

1.006 

.527 

.004 

.74 


The corresponding values of A (obtained from fig. 15 and column (5)), 
T — T 

— r ■ — ■ P (from eq. (All)), and — from the equation 

L r " i p h 


_h 

h' 



(A18) 


(obtained from eqs. (A4) , (A8), and (All)) are also shown in the table. 


The values of in column (j) are actually the first approximation in 

an iterative procedure since h ! and T p * were used to obtain the values 
of A. The final results, of course, depend directly on the accuracy of 
the surface temperatures of the plug and model wall. It is believed that 
the measured values of T v (as shown for this example in fig. 2(d) and 
used in column (5) ) are too large because of the type of thermocouple 
installation, which was described in the section entitled "Hemisphere- 
cylinder heat-transfer model.' 1 The magnitude of this error can be esti- 
mated by assuming one-dimensional heat flow through the model wall of 
thickness Ar m . The relation for the outside-wall temperature T v 
(assuming h = h f ) is then 


Ar 

T. + — -E h'T 
1 k st 


1 + 


Ar 


m 

Sst 


r 


V - 
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which was used to compute the values shown in column (D . These esti- 

mated values were then used to compute a new set of values for — given 

h 1 

T_ 

in column (9) . For this particular test, then, the correction factor — 
varied from about 0.53 to 0 . 78 * 

In view of these possible errors in T w and Tp and the conserva- 
tive effect of the boundary conditions (eqs. (A13)), the correction factors 
shown in column ( 9 ) are believed to be somewhat small. Consideration of 
all pertinent effects indicates that most of the heat-meter data in this 
report should be multiplied by a correction factor of 0.8 with a probable 
uncertainty of ±10 percent. However, when h is small, as at 9 = 90° 
on the hemisphere, the values of Xp are small and larger correction 
factors may be required. It was impractical to attempt to apply a cor- 
rection to all the data because of possible unknown errors in T w> T p> 
and e. 
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(a) Isothermal heat-transfer model. 

Figure 1.- Sketches of models and instrumentation used in the investigation. All dimensions are 

in inches unless otherwise indicated. 
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(b) Equilibria- temperature model. 
Figure 1.- Continued. 
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(c) Hemisphere-cylinder heat-transfer model and heat meter. 


Figure 1.- Concluded 
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(a) Recorded temperature variation with time on the inside and outside surfaces at about the 
same 9 station. ^ = 2.0; p Q = 98 lb/sq in. abs; T = 492° F; = 3.4 X HP (rough 

surface) . 


Figure 2.- Temperature distributions on the heat-transfer models for typical tests. 
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(b) Temperature-time history of the copper sphere as read directly from 
the records. Flagged symbols denote temperatures measured at the 
inside surface; all other values are measured about 0.05 inch from 
the outside surface. ^ = 2.0; p Q = 98 lb/sq in. abs; T Q = 1+92° F; 

= 3*4 x 10 (rough surface). 



(c) Temperature time history of the copper sphere for the same test 
shown in figure 2(a) but corrected for the effect of thermocouple 
depth. 


Figure 2. - Continued. 
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(d.) Temperature distribution on the inside and. outside surfaces of the hemisphere- cylinder heat- 
meter model. = 4.15; P 0 = 324 lb/sq. in. abs; T Q = 241° F; R*, = 3.2 X 10 6 (polished 
surface) . 


Figure 2.- Concluded. 
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Figure 3 .- Distribution of the pressure coefficient around the sphere 
^ = 2.00 for various stream Reynolds numbers. 
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4.- Variation of the separation point with stream Reynolds number on two of the models 

at ^ = 2.0. 
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Figure 5 .- Schlieren photographs of the recovery-temperature model at 
= 2.0 and three stream Reynolds numbers. 
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Figure 6.- Distribution of the pressure coefficient around the sphere 
at = 4.15 for two stream Reynolds numbers. 
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(a) 0 = 0° to 55°. 


Figure 7«- Variation with stream Reynolds number of the ratio of the 
equilibrium temperature to the stagnation temperature at = 4.15 
for two wall thicknesses. Flagged symbols denote negative values 
of 0 . 
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Figure 7*- Continued 
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Figure 7-- Continued 
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Figure 7*- Continued 
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Figure 7« - Concluded. 
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Figure 8.- Variation with 0 of the ratio of the recovery temperature to the stagnation temper- 
ature at = 4.15 for three Reynolds numbers . 
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Figure 9 .- Variation of heat-transfer parameter with wall-temperature 
parameter for various values of 0 on the heat-meter model at 
tio = 4.15 and « 5.5 x 10 6 . 
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Figure 10.- Ratio of the equilibrium temperature to the stagnation temperature at ^ - 2.00 

and several Reynolds numbers. 
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Figure 11.- Variation of Nusselt number at the stagnation point with the product of Reynolds 

number and the dimensionless velocity gradient. 
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Figure 12.- Variation with 0 of the ratio of the local heat-transfer coefficient to the value 
at the stagnation point at = 2.0 and several stream Reynolds numbers. 
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Figure 13 . - Local Nusselt number divided by the square root of local 
Reynolds number at 1^ = 4.15. The solid symbols denote tests made 
without polishing between tests. 
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Figure 15 . - Calculated variation of a correction factor which is used to obtain the ratio of 
corrected to indicated heat-transfer coefficients for the heat meter. 
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